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FORWARD 


The following manuscript is presented in two parts. Part I is 
the discussion of a temperature measurement system which was used in 
the experimental program discussed in Part II. The two discussions 


are presented separately since the temperature measurement system is 
not solely intended for the present experimental program; i.e., it is 
a system that will bo quite useful in many other experimental pro- 
grams. Those details which sre specific to the present experimental 
program will be found in Part II. 


This document is. with one important exception, taken from the 
M.S. thesis of the senior author (C.E. Wark) . At the time of the 
thesis preparation, it was not recognized that an experimental pro- 
cedure error {a constant volume flow rats and tot a constant mass flow 
rate was used from the metering nozzle to the jet exit piano) had 
introduced an error into the computation of J for the heated flow con- 
dition. The J values, of the thesis, were thought to be constant 
(16,64) whereas their true values ere shown in this report. 
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work rate 


streaswise coordinate (see Figure 4) 

transverse coordinate (see Figure 4) 

streaswise lab coordinate (see Figure 11) 

transverse lab coordinate (see Figure 11) 

lab coordinate (see Figuro 11) 

density 

time constant 

kinematic viscosity 

jet coordinate (parallel to Y) (see Figure 11) 
jet coordinate (tangent to jet trajectory) (see Figure 11) 
jet coordinate (perpendicular to § and t)) (see Figure ll) 
non-dimensional temperature (Kamotasi and Greber (1971)) 


mean value 


PART I 

TEMPERATURE MEASUREMENT SYSTEM 


CHAPTER 1 
INTRODUCTION 


1,1 Temperature Measurements 


Accurate temperature measurements have become increasingly impor- 
tant in many areas of combustion, heat transfer, and fluid mechanics 
research. The relevance of such measurements to combustion and heat 
transfer research is quite obvious whereas, in fluid mechanics 
research, one often desires to obtain information regarding a fluid 
flow problem by using heat as a passive contaminant to thermally mark 
the region of interest in the flow. 

Several methods for measuring temperature exist. Bimetallic and 
liquid-in-glass thermometers utilize the familiar process of thermal 
expansion. Optical techniques include Raman Scattering and Rayleigh 
Scattering; both of which require sophisticated equipment (Peterson 
1919 ). Another method uses thermistors which are in the class of sen- 
sors whose resistance changes with temperature. A thermoelectric 
sensor, (thermocouple) is composed of two dissimilar materials which 
form a closed ciruit. If the two junctions sense different tempera- 
tures, an electromotive force (emf) will be established. This effect. 
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which was discovered in 1821 by T.J. Seebeck, is the basis for ther- 
moelectric temperature measurement. In 1834, J.C.A, Peltier 
discovered that when current flows across the junction of two metals 
heat is either released or absorbed at that junction. This can, 
therefore, cause a junction to be at a temperature different from that 
of its surroundings. In 1851 W. Thomson found that an additional emf 
is produced when there is a temperature gradient along a homogeneous 
wire. The total emf generated in a thermocouple is due to both the 
Peltier emf which is proportional to the difference in junction tem- 
peratures and the Thomson emf which is proportional to the difference 
between ihe squares of the junction temperatures. These proportional- 
ity constants depend on tha two dissimilar materials used, and can in 
principal be found from a calibration of the thermocouple with a known 
temperature. In actual practice however, the total voltage is of 
interest and therefore the emf's generated by the Peltier and Thomsen 
effects are not determined separately. It should be noted that irre- 
versible Joule heating effects (I a R) are present and will raise the 
temperature of the circuit above the local surroundings , however if 
the current in the circuit is that due to the thermoelectric effect 
only. Joule heating can be neglected (Doebelin 1 966 ), 

Errors involved in using thermocouples may result from radiation, 
conduction, or time constant effects. These effects cause the meas- 
ured temperature of the probe to be different from that of the 
surrounding medium. The radiation error is a result of the difference 
in temperature caused by radiation from the wire or prongs to the sur- 
rounding walls. Conduction of heat, between the thermocouple wires 
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and the prongs, will introduce error into the temperature measurement 
' This error being known as the conduction correction. The response 

time of the thermocouple can also result in temperature measurement 

J 

errors and, of the above three effects, has received the most atten- 
tion in the available literature. 

Radiation from the probe is proportional to the fourth power of 
the temperature, and thus will have a correction which increases shar- 
ply with temporature. As pointed out by Bradley and Matthews (1968), 
the magnitude of the correction decreases as wire diameter decreases. 
Moffat (19S8) performed a series of tests to determine the effect that 
various levels of radiation loss had on the response time of a 16 gage 
round-wire loop junction thermocouple. For gas temperatures up to 
1600°F it was found that radiation had a small effect on the response 
of the probe. 

The distance (1) from the temperature sensing junction to the 
support prongs is a critical factor if conduction effects are to be 
negligible. Petit et al.(1982) use the ratio 1/ 1 c » 

(21/d w ) ( (kg/k^) Nu)^^ to determine the relative importance of conduc- 
tion effects; specifically, if l/l c > 3.0 then conduction effects are 
assumed to be negligible. Bradley and Matthews (1968) use the ratio 
l/d w to characterize conduction effects. The authors do not give a 
. value for this raftio but cite some specific examples. 

The response time of a thermocouple, also known as its time con- 
stant (A), is a function of both the heat transfer coefficient (h) , 
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and the surface area of the thermocouple bead (A s ), The heat transfer 
coefficient is in turn a function of velocity and temperature. The 
most widely used methods for the determination of the time constant 
are the internal heating techniques. Such techniques as proposed by 
Shepard ead Warshawsky (1952) and Kunugi and Jinno (1959) have been 
employed by Petit et al, (1982), Lockwood and Moneib (1.979) and Bal— 
lantyne and Moss (1977), These techniques determine the time constant 
based on the transient behavior of the thermocouple when the wire is 
run with a periodic current or after the removal of a DC current. 

Petit et al.(1982) made comparisons between the internal heating 

technique of Shepard and Warshawsky (195?) and an external heating 

technique described in Petit et al. (1979). A 3, 5pm dia. Pt— Rh 10*fc 

wire is subjected to • sinusoidal or square current and it located 

just upstream, in a known flow field, from the given thermocouple. 
Their results showed significant differences in the measured time con* 
stants between the internal and external heating techniques, which the 
authors attributed tc both Joule and Peltier effects. 

Compensation for the time response of the thermocouple probe can 
be performed using either analog or digital techniques. The basis for 
which electrical circuits compensate in part for the time response of 
the probe was formulated by Shepard and Warshawsky (1952) . The 
electrical circuit requires a prescribed time constant for use in a 
differentiating circuit. Ballantyne and Moss (1977) studied the sig- 
nal distortion resulting from the use of an average time constant 
value in the electrical circuit. The authors found that where there 
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were large local variations in the time constant values, significant 
distortion of the compensated temperature statistics resulted. In 
regions of moderate heat transfer fluctuations however, the analog 
technique would seem to work well if one could obtain a correct iver 
age value for the time constant. Lockwood and Moneib (1979), in 
measuring temperature fluctuations in a heated round free jet, 
obtained good results u*ing this method. Digital compensation 
requires fast sampling rates and a moderate amount of processing time. 
The main advantage of using digital compensation is that an average 
value for the time constant need not be used, instead variations of X, 
if known, could be implemented into the processing software. 

Of the two junctions in a thermoelectric circuit, one is common^ 
referred to as a reference junction. The reference junction must be 
kept at a known temperature if one is to determine the temperature at 
the desired location. An ice bath is very common and its temperature 
is reproducible to about ,001®C (Doebelin 1966). Another method is io 
allow the reference junction to sense ambient temperature; the net emf 
thus generated would thon correspond to a temperature difference with 
respect to the ambient. 

When using a single thermocouple, one can obtain statistical 
quantities of temperature at discrete locations in space or a tempera- 
ture time series at a given point in space. An array of several 
thermocouples will provide for the above mentioned quantities and, if 
all thermocouples are sampled simultaneously, will also provide for 
instantaneous isothermal maps of the given flow field. The present 
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temperature measurement system is a multi-point, fast response, simul- 
taneously sampled array of 66 thermocouples, which can be used to both 
evaluate the temporal evolution of the instantaneous temperature field 
and obtain temperature statistics at discrete values in the flow 
field. 

1.2 The role of the Time Constant in a Temperature Measurement System 
1.2.1 The First Order Response Equation 


The time constant (1) is the parameter in the first order 
response equation which describes the relationship between the true 
temperature value at the thermocouple location (T ) as a function of 
the measured temperature value (T^) , This equation can be derived 
from a thermal energy balance for the thermocouple bead; namely 

-pcV-3T B m - hA s (T n -T c ) + k a A x ST/dxj^-k^ 3T/8x] Xj (eq 1.1) 

Where V- and A s are, respectively, the volume and surface area of the 
bead. k a and k a are the thermal conductivities of the wires connected 
to the bead at x a and x a , A x is the cross sectional area of the wires 
and h is the convective heat transfer coefficient. 

If the wire is sufficiently long with respect to the diameter of 
the wire (Petit et al 1982 and Bradley and Mathews 1968), then conduc- 
tion is negligible. Equation 1.1 then reduces to 

-po¥idT ffl / dt) - hA s (T m -T 0 ) (eq. 1.2) 
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■which can be written in the form 


X(dT n /dt) + T a “ T c (eq. 1.3) 

where X“peV/hA |4 For a givon wire material, i.a. for given values of 
p and c, it is seen that X depends upon the ratio of the bead's volume 
to surface area (¥/ ) , and upon the reciprocal of the convection 
coefficient, h“^. From convective heat transfer considerations, h is 
a function of the velocity and the thermal properties of the fluid. 
The latter may be expressed as a function of temperature for the given 
test section pressure, which is essentially constant for the small 
overheat values of the present study. Considering the former, h may 
be approximated as being proportional to u ' for a laminar boundary 
layor; hence, the time constant X would be expected tc be proportional 
to 

1.2.2 The Hole of X in the Data Acquisition Process 


The true temperature (T.) may be determined from the measured 
temperature (T m ) if the latter measurement is supplemented with a 
simultaneous evaluation of X and dT ffl /dt. Specifically, consider that 
a high-rate time series of T m (t) values is available and that the time 
derivative (dT m /dt) is evaluated from these data. T c can then be cal- 
culated from (1.3). The data samples, so provided, could be used to 
form ■ sample population of the temperature values that exist in the 
flow field. The statistical values of interest can then be determined 
from this sample population. 
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The genetil data acquisition procedure, that is considered 
herein, is for three samples to be acquired in rapid sequence: 
^ra^k-l' ^m^k' T m^k+1' *Hw* dT ffl /dt to be defined from a central 

difference scheme and the corresponding value of T c )^ to be computed. 
The T c )j. values are retained as the primitive data for further evalua- 
tion. 


It is noteworthy that the mean value, of the T c values 30 
defined, would be independent from- X under a relatively 

non-restrictive assumption. Specifically, consider that the velocity 
magnitude (u) is simultaneously measured with the T ffl value. The cor- 
rected temperature (T c ) is then given by 1,3 where X=Mu) is used for 
the time constant. Using an over-bar to denote the mean of the sample 
population and a prime to denote the deviation (i.e,, difference 
between a given realisation and its mean value), the eqnation for the 
mean value of the corrected temperature of the sample population is: 

T c = \ + XUTjJtT + IuT m /dt5 (eq. T.4) 

For 1 stochastic process, (dT m /dt)*0 and T c =T n if there is no correla- 
tion between the X and (dT m /dt) deviations. 

The variance of the sample population is, however, clearly 
affected by the time constant value. This can be shown as follows. 
Subtracting 1.4 from 1.3 yields 

T c>k 53 Vk~ \ = Vk“ *n + X k (dT m /dt) k" i<4T B /dt} (eq 1.5) 


.55 
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It vill be assumed that the correlation: i(dT m /dt) is zero for 
this analysis, Given this assumption, the variance of tho sample 
population is 

2'^) k X k (dT 01 /dt k )+^(dT n /dt)* (eq 1.6) 

The second term, r.h.s. of (1.6) is not expected to have a signifi- 
cant value; the third term is, however, always positive and hence the 
sample variance for the corrected temperature will exceed that of the 
measured temperature if X > 0. 

The following chapters of Part I identify the physical tempera- 
ture measurement system and the data acquisition scheme. Chapter 2 
provides results for, and a description of the method used to deter- 
mine, the time constant values for all thermocouples. A brief 
discussion of the velocity dependence of the time response is also 
included. The measured temperature value is digitally corrected for 
the time constant of the thermocouple; this correction scheme is 
explained in Chapter 3. 


CHAPTER 2 




TEMPERATURE MEASUREMENT SYSTEM 


2 .1 THE PHYSICAL SYSTEM 


2.1.1 Thermocouples end Support Array 


A schematic representation of the present thermocouple probe is 
shown in Figure la. The diameter of the probe body is 1.58mm and the 
overall length of the probe is 153mm. The micro-fine wire thermocou- 
ples (7.62 |im dia.) are supported by ,235mm dia. prongs separated by 
1.59mm (see Figure lb). The thermocouples used in this study were 
type E which is made up of a positive conductor of nickel-chromium 
(chromel) and a negative conductor of copper-nickel ( Constantsa) , the 
sensitivity of the type E thermocouple is 42.5 ji volts/*F (Doobelin 
1 966 ) . 

The 72 thermocouples are supported in the array shown in Figures 
2 and 3, The frame is 145mm square. There are 18 modular elements 
with four thermocouples supported in each modular element as shown in 
Figure 3. The 18 modular elements are connected by ''forks" with two 
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of the corner element* secured to the frame. The area blockage aS the 
measurement point is 1.5%, and at a downstream distance of 102mm, the 
area blockage is 42%. Figure 2,14 in Schlicting (1968) shows that the 
affect on the velocity field at the measurement point would be negli- 
gible even if the forks and modular elements created a 100% blockage. 


2 . 1.2 Reference Junctions and Assembly 


The 72 reference junctions for the thermocouples are located in 
the boundary layer bleed passage of the Free Shear Flow Laboratory 
(FSFL) flow facility (Figure 4). The reference junctions sense T^ an d 
they have a relatively large time constant since the wire diameter is 
large (,8mm). The thermocouples and reference junctions are connected 
via 72 twisted pair shielded cables. The cables are 35 feet in 
length. Extra lengths of chromel and constantan were welded on to tho 
leads of the thermocouple probes, as shown in Figure 2 , to ensure 
that no second thermocouple effect would be present in the measure- 
ments, i.e., if the chromel/ copper or constantan/copper junctions were 
subjected to a temperature above the ambient value, T,,, an additional 
emf would be induced. 


2.1.3 Electronic Processing Units 

A schematic reprsentation of one of the four electronic process- 
ing boards is shown in Figure 5. Each processing board has 20 sample 
and hold units (detail shown in Figure 6), one multiplexer with 32 
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channels and 20 amplifiers. The resistors used to amplify the thermo- 
couple signal were metal film with a 1% tolerance; their values were 
selected to give a nominal gain of 500. The reference junctions were 
mounted on the back of the processing boards thereby reducing 60 Hz 
noise pickup by the leads connecting the reference junctions and pro- 
cessing boards. Placement of the processing boards in the FSFL flow 
facility can be seen in Figure 4. 


The processing boards were designed and built by Sigmon Electron- 
ics. 


2,1.4 Oats Acquisition System 


The input to the processing boards is a voltage corresponding to 
the temperature difference between the thermocouple probe and the 
associated reference junction, T m -T m . The signal is amplified, sam- 
pled and held, multiplexed and presented by the processing boards to a 
four channel TSI 1075X. analog to digital converter. The multiplexers 
on the processing boards are controlled by a DEC DRVll interface board 
which allows the user to pick the starting channel and the number of 
channels on each board to be read. The 12 bit analog to digital con- 
verter (A/D) has a 0-+5 volt range; therefore, the resolution of the 
A/D is 1.2 millivolts. By connecting all channels of the A/D to a 
common signal, it was found that the accuracy of the A/D was ±2.4 mil- 
livolts. The four channels have matched 50kHz input filters and are 
sampled simultaneously. The A/D is controlled by a direct memory 


access interface board (DRV11B) which is installed in a Charles River 


Data Systems PDP1123 microcomputer . 
2.2 System Performance Evaluation 


2.2.1 Evaluation of the processing boards 


Several pertinent characteristics of the processing boards were 
evaluated prior to the connection of the thermocouples. The charac- 
teristics checkod were i) the gain for each channel, ii) simultaneous 
sampling of all channels and iii) the sampling rate threshold of the 
multiplexers. Since all 80 pairs of resistors are not identical, the 
amplification of the thermocouple voltages will not be identical for 
all channels. The evaluation of the gain (G) for each channel was 
performed as follows. Known d.c. voltages (E^) were applied to the 
inputs (points 2,3 of Figure 6 ) for all channels and the previously 
discussed data acquisition system was used to sample the output (E Q ) 
of each channel. For the transfer function; 


E o - G ‘ E i - E „ff> 


(eq. 1.7) 


where E off = offset voltage, the gain (G) was obtained by a linear 
regression fit to Equation 1,7 for the data p«.irs: (E^,E 0 ), The 

offset value, determined by a separate procedure as described in 

chapter 3 ; was performed by giving the channels a sequence of DC vol- 
tages and using the previously discussed data acquisition system to 
sample all channels. A linear relationship exists between the thermo- 
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couple reading and the input DC voltage; namely, the thermocouple 
reading is equal to the channel offset plus the DC voltage multiplied 
by the channel gain. By using a linear regression routine, the value 
for the gain for every channel was found. 

To determine if all channels were sampled simultaneously, a com- 
mon signal was input to each of the channels and the channel outputs 
were compared. Specifically, a nominally 200 Hz sine wave from a sig- 
nal generator was applied to the channels and a visual comparison of 
the outputs, from various channels, showed that the sample and hold 
units were sampled simultaneously. 

It was observed that the multiplexers would sometimes get "lost" 
at high sampling rates; i.e., it seemed that the multiplexers on the 
processing boards did not start on the user specified starting chan- 
nel. This hardware problem was compensated for by the data 
acquisition software, as described in the following chapter. 


2.2.2 Thermocouple Time Constant Evaluation 


The experimental configuration, to determine the time constant 
values, consisted of placing a 5pm diameter tungsten hot-wire probe 
directly upstream of the subject thermocouple. The hot-wire probe was 
positioned easily and precisely by means of a traversing mechanism. 
This is shown schematically in Figure 7. 
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The data acquisition system discussed above was osed to sample 
the subject thermocouple and the hot-wire output. During the time 
period that the A/D is sampling, the hot-wire anemometer is switched 
from operate to stand-by, thereby causing sn impulsive decrease of the 
heating current to the hot-wire. The hot-wire has a sufficiently fast 
response that the thermocouple essentially senses a step change in 
temperature. Approximately five runs at a constant flow speed (~ 
7mps) were performed for each thermocouple. 

A schematic representation of a thermocouple's respond (with 1, - 
5msec) to the hot-wire can be seen in Figure 8. The initial flat por- 
tion of the thermocouples' response is attributed to the separation 
distance between the hot-wire and thermocouple, and the convection 
velocity. The asymptotic value of the thermocouples' response would 
be the laboratory temperature (T^) for this case. 

The solution to equation 3 can be written as follows: 

ln[(T m -T.)/(T r T^)] - -X -1 t (eq. 1.8) 

The terms in equation 4 were evaluated from the time constant data as 
follows. T^ was the initial value of the thermocouple just before the 
hot-wire ■nemometer was switched from operate to standby, T ffl was the 
measured value at the sequence of sample times t, and T^ was the 
steady state valuo reached by the thermocouple after the passage of 
the trailing edge of the hot-wires' thermal wake. A linear regression 
routine was then used to determine X <* -slope - *) when the data were 
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fit to equaton 1.8 . 


An average value for A was determ. * ;d based on nominally five 
data rnns per thermocouple. Consistent values of A ,per thermocouple, 
were found when the hot-wire probe was placed directly upstream of the 
thermocouple probe. This proper placement resulted in a relatively 
large T <) -t w (‘*24 °F) and a low standard deviation of the data when fit 
to equation 4. Oscilloscope traces of the thermocouples response were 
also observed to be very smooth. 

The procedure described above proved to be a good check on the 
integrity of the thermocouples. This check revealed 66 good thermo- 
couples, for which the variability of A at 7mps was found to be 2.0 
msec to 4.7 msec. A nominal uncertainity of -to. 20 msec is associated 
with the A value for a given thermocouple. 


2.2.3 Reynolds Number Dependence of Time constants 


In addition to the data runs at 7mps, time constant values for 
four thermocouples were also found for the range of velocities: 2mps 

to lOmps. The A(Re) results can be found in Figures 9 and 10. Figure 
9 is a representation of A vs. u and Figure 10 is a plot of A vs 
u~l/2 ^ere u is proportional to Re for a given diameter and the nomi- 
nally constant kinematic viscosity values (6% difference between v_. w 

max 

tnd v min) • 
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The X(u) results appear to approach an asymptotic value of X for 
large u. The experimental results of Figure 10 support the observa- 
tions made in 2,2.2; specifically that X is proportional to u" 1 / 2 . 


CHAPTER 3 


DATA ACQUISITION PROCEDURE 


3.1 Offset De termination 


The voltage, that is sampled by the A/D represents the net emf: 
re f ] , pins the offset associated with the electronic com- 
ponents for that channel. The offset for each channel most be 
subtracted from the voltages sampled by the A/D before any other pro- 
cessing steps can be performed; therefore, determination of the 66 
offsets is an important part of the data acquisition procedure. 

By placing the thermocouple array in the ambient temperature flow 
field of the FSFL flow facility, the temperature sensed by the thermo- 
couples (T m ) and reference junctions (T n ) will be approximately equal. 
Discrepancies between the two temperatures will exist insofar as there 
are thermal variations within the room; these thermal variations will 
hereafter be referred to as “thermal noise". The voltages sampled by 
the A/D will then simply be the offset values for the 66 channels, 
note that [E t c (y,z,)- E re £]=0 for offset determination. It is 

important to note that a stabilization period is required for the 
offsets. Apparently this stabilization period is needed to bring the 
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processing boards to a ataady state temperature value once the boun- 
dary layer bleed auction is activated, 

3.2 Sampling Rate Considerations 


Referring back to equation 3 one can see that the quantity dT m /dt 
it essential in determining the corrected temperature value (T c ). 
Evaluation of dT m /dt is performed using central difference techni- 
que, therefore a high sampling rate is required. 

When choosiRg a sampling rate, the rate threshold of the 
multiplexers must be considered (refer to Part 1 Chapter 2). It was 
found that at the highest sampling rate possible for the A/D(3.125 kHz 
for 80 channels) the multiplexers would get "lost"; i.e. they would 
not start on the user prescribed starting channel and hence, the order 
in which the data was stored in memory would not be correct. At the 
next highest sampling rate (1.56 kHz) the multiplexers would start on 
the prescribed starting channel only a portion of the time. Therefore 
a check is performed in software to determine if the order that the 
data is stored in memory is correct. If correct, then the data run 
will be further processed (as will be explained in the following sec- 
tion) and if incorrect then the values are deleted and another data 
run is taken. The "lost" condition of the multiplexers improves as 
the sampling rate decreases and the user can choose one of six sam- 
pling rates between 20 Hz and 1.56 kHz. 

3,3 Time Constant Correction and Data Storage 
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Two modes of data storage are available to the user of the 
present temperature measurement system, A time series of data is 
required for studying the instantaneous behavior (e.g. instantaneous 
isothermal maps) of the flow field in question, and a histogram, com- 
posed of statistically independent samples, is used to determine 
statistical quantities at the discrete locations in space occupied by 
the thermocouples. 

The data which are sampled and stored as v j time series, are taken 
at a high sampling rate (1.56kHz) so that the measured temperature 
value can be corrected as in equation 1.3. This high sampling rate 
(,64msec) is appropriate when considering the time constants of the 
thermocouples (=4 .5msec at 2mps) The value for the time constant used 
in equation 1.3 is baaed on a user defined velocity or it can be based 
on a user defined function of velocity with temperature. As explained 
in the previous chapter, X values were found for all thermocouples at 
a convection velocity of 7mps. Using these results, and if one 
assumes a common slope of - 5 (msec) / (mps)"*-' (Figure 10) for all ther- 
mocouples, the value of X for any velocity in the range 2 1 u i lOmps 
can be found for each thermocouple. 

As noted in the Introduction, the first order correction of the 
measured temperature values needs to be performed in order to obtain 
accurate values for statistical quantities of the data. The corrected 
temperature value T c , is used to build a histogram for each thermocou- 
ple. The program is written such that the user defines the minimum, 

the maximum and the number of intervals for the T„ distribution. For 
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each data run a total of tavea points pa t thermocouple are sorted Into 
the correct interval of the approprir.ve histogram. Although thore are 
approximately 96 kBytes of memory used for the measured temperature 
values per data run, only seven statistically independent samples are 
kept. This seemingly excess amount of raw data is required because of 
the high sampling rate (.64 msec) necessary for the central difference 
scheme. The user specifies the number of points desired to develop 
the histograms and the appropriate number of data runs will be taken 
under software control. 


PART II 

JET IN A CROSS FLO? EXPERIMENT 
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CHAPTER 1 
INTRODOCTION 


A jet exiting into a cross flow is an example of a complex three 
dimensional flow field. The interaction of the jet fluid with cross 
stream fluid results in two large counter rotating structures which 
are one of the dominant features of this flow field (Figure 12). 
These structures are responsible for the distribution of entrained 
cross stream fluid within the jet. 

The interaction of the jet with cross stream fluid results in 
high pressures on the upstream surface of the jet column. The down- 
stream face of the jet is characterized by lew pressure values and 
strong entrainment of the ambient fluid. A qualitative discussion of 
the flow field, including surface topology effects, is presented by 
FossURRO) . 

The deflection of the jet makes it useful to define a system of 
jet coordinates 5 follows the jet trajectory, n is perpen- 
dicular to £ and parallel to Y, and s' is perpendicular to both £ and 
n* Figure 12 shows this system of coordinates and the lab coordinates 
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X, Y, and Z. 


The important parameter in this problem is the momentum flux 
ratio (J»p jVj*/p a) U 0 *) . When the cross stream and jet fluid have the 
same density, the velocity ratio (R-Vj/U 0 ) is often used to character- 
ise the problem. 

Some of the physical applications of a jet exiting into a cross 
flow are i) the discharge from a chimney stack, ii) waste discharge in 
rivers, and iii) the cooling of combustion gases. The U tter is the 
motivating application for the present study. 

The study of a jet exiting into a cross flow has been the subject 
of numerous investigations. Keffer and Baines (1962) present velocity 
measurements at a large number of points in the flow field, for the 
velocity ratio values Vj/U 0 of 2, 4, 6, 8, and 10. Isobars for the velo- 
city ratio values of 2.2 to 10 can be found in Fearn and Weston 
(1975). Ramsey and Goldstein (1971) used a large diameter thermocou- 
ple (,13cm) and obtained temperature profiles of a heated jet exiting 
into a cross stream at low PjVj/pmCj,, values (.1 to 2) 

Kamotani and Greber (1971) provide a substantial foundation for 
the present work. Their study was also motivated by the combustor 
cooling problem and they provide velocity and mean temperature infor- 
mation in planes that are perpendicular to the jet axis. Extensive 
reference tc the detailed considerations in their report are made 
herein. Androopoulos (1983) in a recent publication, provides infor- 
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mation regarding the scan and fluctuating thermal field for Ji4.0. 

These previous studies have been performed with a relatively low 
disturbance cross flow. The flow field in a combustor is. however 
disturbed and it is this cross stream disturbance effect on the dynam- 
ics of a jet in a cross flow which is the primary focus of the present 
study. 

The present study uses thermal energy to mark the jet fluid. The 
temperature measurement system described in Part I is used to investi- 
gate the temperature field downstream of the jet exit. The cross 
stream disturbance effect is investigated by comparing the results 
obtained for a jet exiting into a disturbed cross flow with those for 
a jet exiting into an undisturbed cross flow for c given momentum flux 
ratio. Two groups of momentum flux ratios: J=17,19,21 and 1=67,76, 
and 84 and three jet overheat values: (T,-T ot> )“22.2,41,7, and 61.1 °C 
were used in the present experiment. 


CHAPTER 2 


EXPERIMENTAL EQUIPMENT AND PROCEDURE 


2.1 Experiment!? Equipment 


2.1.1 Facility for the Crosa Flow 


A schematic representation of the FSFL wind tunnel can be seen in 
Figure 12. This flow facility provides a large.planar shear layer for 
use in the disturbed cross flow case of the present study. 
Specifically the jet exhausts into the free shear layer at a location 
of tt /U re £ = ,25 at the end of the three meter test section. The 
u/U ref and turbulence intensity profiles based on both the free stream 
and local velocities can be seen in Figures 13,14 and 15. The Q /^ re f 
= .25 (y=80cm) location in the free shear layer was chosen based on 
the relatively small velocity gradient and the relatively high turbu- 
lence intensity valv.es existing at that point (Figures 13 and 15) . It 
is pertinent to note that, based on the energy spectrum results at 
y**80cm (Figure 16), the high turbulence intensity values are caused by 
very low frequency fluctuations. For the undisturbed cross flow case, 
a "false wall" added to the FSFL flow facility (shown as a dotted line 


in Figure 4) resulted in « quiescent flow field at the jet exit loca- 
tion. The turbulent intensity at this location was -.6% and the 
energy spectrum results can be seen in Figure 17. 

A boundary layer plate* the leading edge of which has a parabolic 
profile, was positioned 203.2mm below the t«ip of the Test Section 
Extension (T.S.E.) (Figure 18). An adjustable diffuser flap, for the 
parallel alignment of the leading edge stream line, was attached to 
the downstream edge of the boundary layer plate, (Figure 19). 


2.1.2 Jet Delivery System and Exit Conditions 


A Spencer blower is used to provide the necessary flow rate for 
the i »t. The exhaust from the blower is ducted into a 136cm by 101cm 
plenum. This is shown schematically in Figure 18. Flow enters the 
jet. delivery tube through a 6.35mm dia. sharp edge orifice. At a 
downstream distance of 10mm from the entrance orifice, a static pres- 
sure tap is located in the wall of a 25mm I.D. copper pipe. The 
static pressure from this tap, minus that of the plenum, was used to 
measure the flow rate. After the right angle bend, a 50mm I.D, 
galvanized pipe is used to deliver the flow to the jet exit. 

The jet exits through a lOt .1 dia. sharp edge orifice (Figure 
18). Care was taken to ensure that the exit plane of the jet was 
flush with the lower surface of the boundary layer plate. The center 
of the exit orifice was located 260mm downstream from the leading edge 
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of the boundary layer plate, as shown in Figure 19. Figures 20 and 21 
are plots of the cross stream boundary layers, just upstream from the 
jet exit orifice, for the undisturbed and disturbed cross flow condi- 
tions respectively. Also shown in figures 19 and 20, are the 
reference curves for a turbulent boundary layer in a zero pressure 
gradient flow, viz., u/0*, = (y/6)*^, and the Blasius solution for a 
laminar boundary layer. 


2.1.3 Jet Heating Considerations 


A Superior Electric powerstat was used to supply a controlled 
voltage level, to a flexible resistance strip heater. The strip 
heater was wrapped around the copper delivery pipe in the plenum above 
the T. S.E. . Fiberglass insulation, 220am in dia. , was placed around 
both the strip heater/copper pipe and the galvanized delivery pipe to 
protect from heat loss to the surrounding plenum and/or cross flow 
(Fignre 18). A reference thermocouple (wire dia of .254mm) is located 
lOmrn upstream from the jet exit and was used to monitor the jet tem- 
perature. The leads of the thermocouple extend through the sides of 
the pipe and are connected via one of the shielded cables, to a chan- 
nel on the processing boards. 


2.2 Experimental Procedure 


2.2.1 Veloaity end Jet Temperature Determination 


The cross stream velocity, ( ) , was determined, for the 
undisturbed cross flow case, by using a pitot-static pressure probe 
placed in the free stream of the FSFL flow facility* For the dis n ' 
turbed cross flow study, the 0 o value was 1/4 of the velocity 
calculated from the pressure probe reading. 

The mass flow rate through the delivery tube was determined using 
the static pressure (downstream from the entrance orifice) and the 
plenum pressure. The discharge coefficient, (Cj), which is required 
to correct for the non-uniform velocity profile at the orifice, was 
known a-priori as a function of prossure drop across the orifice, A 
thermometer located in the plenum, provided for the density value used 
in the flow rate calculation. The average jet velocity (Vj) was 
determined from the calculated mass flow rate. 

The output from the reference thermocouple, used to monitor the 
jet temperature, is sampled using the previously discussed data acqui- 
sition system. Knowing the sensitivity of the type E thermocouple, 
76.5pmvolts/ # C, an output voltage corresponding to the desired jet 
temperature can be calculated. The desired output from the reference 
thermocouple is obtained by adjusting the voltage level of the power— 
stat. 
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It was found that a period of approximately one hour was required 
for the jet temperature to stabilize, however a period of three hours 
was required for the jet mass flow rate to stabilize. The latter was 
a direct result of the initial transient period of the Spencer blower. 


2,2.3 The Temperature Meaaurement System 


The results of Kamotani and Greber (1971) were used to determine 
an appropriate position of the thermocouple array with respect to the 
jet exit location. Their results include the Vj/u 0 values of interest 
(4 and 8) for the present study. The thermocouple array was posi- 
tioned in X/d (*0.25) corresponding to the chosen £/d value. The Z/ d 
position (ho, 33) was chosen such that the maximum amount of tempera- 
ture field information would be obtained, (The Y/d (*0.25) placement 
was to simply center the array with respect to the jet exit.) The 
array was mounted on a pair of dove tail traverses, allowing for X and 
Y positioning, which were mounted on an aluminum plate positioned in Z 
by a scissors jack (see Figure 19). The measurement plane of the 
present study was located at X/D ** 3.9 for the lower momentum flux 
ratios: JX21 and then moved to X/D = 3 .5 for the higher momentum flux 
ratios: 1161 where the displacement: AX/D, was accurately 
controlled: *0.05. 

The data acquisition procedure, discussed in Part I Chapter 3, is 
used to acquire the data lor the present study. The following 
features are specific to the present experiment, i) The offsets asso- 
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ciated with the channels <s2 the processing board* require a 
stabilization period (aa noted in Part 1 Chapter 3), which was found 
to be approxiaately 43 aiautes for the pretest study, ii) The sta- 
pling .rate of 1.56 kBz <-64asec) is appropriate when considering the 
central difference scheae used to evaluate the tern dT a /dt; that is. 
.123 jet disasters of fluid will be eonvected downstream every 
,64nsec. This estinate of the eonvected distance is bated on a con- 
vection velocity of 2nps which it the vslue of the free streaa 
velocity for all caaea. 


CHAPTER 3 


DISCUSSION OF EXPERIMENTAL RESULTS 


3.1 Introduction to Discussion of Results 


The primary motivation for this experimental program was to 
determine the effects of cross stream disturbance on the jet in a 
cross flow problem. A secondary body of results is also presented; 
the influence of the momentum flux ratio in the ranges: nil ill and 
67ijiS4, on the location of the jet column. This secondary body of 
results will be discussed after the general features of the measure- 
ments nre identified. The influence of the disturbed cross stream 
will then be considered. 

The influence of assuming f.hat time constant values for each 
thermocouple could be taken as X*X(2mps). was evaluated by processing 
one (high-rate time series) data set with the X values: A«A(5mps) and 
A=0. The former represents an upper bound for the expected velocities 
at the X/d plans of the measurements (Kamotani and Greber (1971)) and 
the latter represents the absence of a time constant correction. From 
the considerations in article I. 1.2. 2, it is expected that negligible 
differences in the mean values of the measured (X a 0) and the 
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Alto, it it 


corrected: A(2mp«) tod X(5mp») should bo obtiined. 

expected that the rat temperature value* should vary at T' 1 0.-0) < T 
'* (X(5mps) ) <T'* (X(2mps) ) . 

Dtiag one of the tiae aeries data sett for: J-84, Tj-T.-dl.l^C , 

and a diiturbed cross flow, the results shown in the following table 
were obtained. 


T r ()c) - TgUUmps)) Tr* - Tr* (A(2mps ) ) 



min 

median 

max 

min 

median 

max 

X-0 

-0.63 

0.05 

0.73 

1.65 

6.0 

17.86 

X-5 

-0.15 

0.03 

0.19 

0.84 

2.0 

7.0 

TABLE 

II. 1 

Influence 

of X on measured 

temperature 

values 



Note: The "min", "median", and "max" refer to a given 
thermocouple (l-of-64) in the array. 

The relatively small differences between the 2 and 5 mps tiae constant 
conditions support the inf err once that the use of a constant X value 
does not exert a significant influence on the present results. 

The temperature field results are most oasily interpreted if the 
conduction of energy through the nozzle plate, and the corresponding 
thermal contamination of the boundary layer fluid, can be neglected. 
The magnitude of this contamination can be conservatively estimated by 
assuming that a thermal boundary layer began at the juncture of the 
nozzle plate and the boundary layer plate (L t ) , A measure of the con- 
tamination is provided by the assumed two-dimensional boundary layer 
at the location of the jet centerline. Using the apapropriate heat 
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transfer relationship for this laminar boundary layer (e.g., Gebhart 
(1971)) : 

q * .664 k (l f -T a ) (U„/v) l/2 (L, 1/2 - l^ 1/2 ) (eq. II. 1) 

and equating this to the flux of thermal energy carried by the boun- 
dary layer: 


q « j pc p <T - T„) u dy (oq II. 2) 

one can solve for the "non-dimensional thermal contamination measure": 

J [ (T-TJ/(T W -T.) ] (u/U„) d(y/d) - 

.664 {t/pc v )/W„d)Wjv) 1/2 (<L ,) 1/2 -(L l ) 1/2 ) - 

1.4 x 10"* (eq. II. 3) 

for U e » 2mps and L, -2 60mm and L 4 *> 213mm. Hence, it is inferred that 
the boundary layer heating effect is of negligible importance for this 
flow field. 
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3.2 The rant 1 Energy Considerations 


Applying the energy eqoation to the control volume shown in fig- 
ure 21 one obtains 


6-^ s **d/dt c> y> J(Vj*/2+gz+u)pdV+ c> # ^ J(Vj*/2+gz+u+p/p)pV.ndA (XI. 4) 

Assuming i) steady state, ii) heat-transfer rste and work rate between 
system and surroundings * 0, and iii) that the kinetic energy term is 
small compared with the enthalpy term, (h«p/p+u), equation II. 4 
reduces to 



pCpT V • n dA 



pCpT V‘n dA 


(oq. II. 5) 


The above integral quantities if evaluated, would provide a good check 
on the experimental data; however, velocity information is not avail- 
able from the present study. The present results provide contours of 
constant temporature (isotherms) for both tho mean flow field and the 
instantaneous flow field. These results show that, at the X-location 
of the thermocouple array, both the mean and instantaneous tempera- 
tures existing in thi3 flow field axe substantially lower than the 
initial jet temperature. Namely, the peak temperatures do not exceed 
the value; 0.25(Tj-T„) , Equivalently, this result shows that substan- 
tial diffusion (i.e. molecular transfer) of thermal energy has 
occured between the jet exit and the observation plane since this is 
the sole agent that can reduce the energy stats of a given fluid ele- 
ment. From (II. 5), and for jet velocities which are approximately 
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equal to the cross stress velocity, it is apparent that the contamina- 
tion field at c.s. 2 oust be much larger than the area of the jet 
orifice. This characteristic of the isotherms is quite apparent as 
noted below. 

3.3 Jet temperature considerations 


For the present experiment, it was desired to use thermal energy 
as a passive scalar contaminant with which to mark the jet fluid. The 
passive contaminant assumption is given strong support by the Kaaotoni 
and Greber (1971) study wherein Tj-T^MI . 7°C and 177.8 ®C overheat 
values resulted in a quite small (-2 percent) change in the z position 
of the maximum temperature. 

This small change was observed at the x/d location of the present 
study and the change was opposite to that expected if a buoyancy 
effect were reaponsible. The effect of buoyancy is expected to be 
negligible since the Richardson number (R.« g (T;-T*)d/U?T«.) for the 

* ® J * ^ I 

present study is small (R^O.005), see Townsend (1975). Figures 22-33 
represent the mean isotherms for the momentum flux ratios: (J) of 
lliJill and 57^J^84, and for three different overheat conditions: 
Tj-T <B **22.2,41. 7. and 61.1 ®C for both the undisturbed and disturbed 
cross flow cases. These contours were obtained using the Surface II 
contouring package on the Cyber 750 at Michigan State University. The 
+ signs on the figures represent the positions occupied by the thermo- 
couples. The coordinates of all contours are those of the 
thermocouple array (s and t). (It is pertinent to note that the posi- 
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tion of Che thermocouple array *at not changed as the jet temperature 
was changed £or a given J value . ) Listed in the caption of each figure 
are the transformation equations between the array coordinates and tie 
jet exit coordinates. The numerical values of the isotherms represent 
1000x(T c -T w) / (Tj-T„) - T r . 

Figure 34 is a plot of the s component of the centroid as a func- 
tion of momentum flux ratio for the several cases investigated. The 
centroid represents an integral measure of the temperature field; in 
this sense it represents s smoothing of the discrete data samples. 
The centroid waa calculated by the equation 

a c.g.“} *I A I / } A I * t c.g.“^ fc I A I 1 1 A I (,q * II * 6> 

where Sj is the s component of the centroid for each isotherm and Aj 
represents the area of each isotherm and was calculated using the 
digitizer on the Prime computer at Michigan State University. 

The results, shown in Figure 34 are quite instructive. Firstly, 
the local variation in J, created by varying T- with a constant volume 
flow rate into tho heater tube, shows that the resulting 
centroid-momentum flux relationship is consistent with the global var- 
iations created by the increase in the volume flow rats. Secondly, a 
consistent pattern, from the four (experimentally independent) sets of 
data, is revealed wherein the jet penetration for the disturbed cross 
stream is less than that for the undisturbed case at the same J value. 
(Incidentally, the consistency of this pattern suggests that the 


uncertainty estimate (±0.7) Cor the s-z transformation is too conser- 
vative.) As noted on che figure, the X/d locations of the two groups 
of momentum flux values are slightly different; it is inferred that 
the trend of increasing centroid distance, with respect to the J 
vslue, would be even more steep if the correction for the 11% change 
in X/d could be accounted for in this graphical representation. The 
influence of the isotherms, selected for the centroid evaluation, ii 
seen to be a second order factor in the above interpretations. (The 
lack of a complete set of isotherms, within the space of the thermo- 
couple array, required the calculation of the two "centroid" values 
that are presented in Figure 34 for the large J values.) 


An alternative method of calculating the centroid is to consider 
the measured temperature as a discrete "mass" point and to compute the 
centroid of the distributed field of such "mass" points. 
Specifically, the centroid was calculated using the equation 


s 
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(eq. II. 7) 


where s^, are the coordinates for a given thermocouple. The results 
cl this technique are not considered to be as valid a measure of the 
thermal field; however, this calculation is dramatically easier to 
execute and it will be employed where the use of the digitized isoth- 
erm contours would be impractical; see Section 3.5. The centroid 
results, using this method, are not presented herein. However, it is 
pertinent to note that the centroids, inferred from the isotherms. 


show a much larger variation, and a more regular dependence upon, the 
momentum flux values than the centroids calculated from the discrete 
mass points. 

A noteworthy aspect of the mean isotherms is the trend: "the 
maximum mean temperature decreases with increasing overheat". This 
decrease in maximum mean temperature can be seen in the mean isotherms 
for both the disturbed and undisturbed cross flow conditions (Figures 
22-33). 

The momentum flux ratio influence can also be seen when looking 
at the histograms obtained for all thermocouples in the thermocouple 
array. A selected few of these histograms can be seen in figures 
35-42. Figures 35. 37, 39, and 41 represent the histograms for a thermo- 
couple located in approximately the center of the jet while figures 
36,38 ,40, and 42 are the histograms for thermocouples located farther 
from the boundary layer plate. 

3.4 Comparison of Present Results with Kamotani and Greber 


Significant differences can be found when comparing the mean 
isotherms of Kamotani and Greber (1971) with those of the present 
study. The absence of kidney shaped isotherms is one factor albeit 
this may relate to the vertical vs, the j et-axi 3-normal plane of the 
data. Compounding this difference in observation planes is the effect 
of an increased jet penetration and therefore the steep intersection 
between the jet column and the X/D“ constant mearsnr ement plane. The 


"steep intersection" is described in the following. 


In order to compare the present results with these of Kamotani 
and Greber, it is required to provide a common momentum flux ratio as 
well as the stme 1/D location. Their comprehensive data for Z m ( x ) 
(i.e,, the Z location of the maximum temperature) readily allows the 
latter; more complicated steps ajje required to provide the "same" 
momentum flux ratio. The apparently linear variation of Z Q ^ with 
respect to J (see Figure 34) allows an approximate interpolation 
between the measured conditions, of the present study, to obtain a 

Z cg _ value that would correspond to the R=8=Vj/Up value of the Kamo- 
tani and Greber investigation. 

Their equivalent J value is inferred to be; 

J = PjVj* / PpU,* 

= (T 9 / Tj )R i 
= ( 530/ 600)8 3 - 56.5 

and the corresponding Z c ^ value (from Figure 34) of the present 
study is 12.5. 

The minor effect of the different X/D values (3.R cf 3.5), was 
neglected in arriving at this value. A ratio of 1.03 was noted for 
( Z m/Z c< g c ) for the present data. Combining these observations allows 
an estimated Z m value to be added to the Kamotani and Greber center- 
line distribution; see Figure 43a. It is noted that the present jet 
penetration is much greater than for their study. A similar procedure 
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coaid not bo readily executed for the lower R value: 4, of their 
study since it would have required an extrapolation from R=19 to 
K*14,l on Figure 34. However, the inverse process: adding their data 
to a nodi f led Figure 34. could be executed and the results are 
presented in Figure 43b. For this calculation, the ratio was 
1.11. It is again noted that the Z c>8> values are substantially 
greater for the present investigation. 

It is inferred that the boundary conditions are responsible for 
depths for the strong differences in the penetration depths for the 
two studies. Specifically, the increased penetration of the jet for 
the present study is attributed to the higher effective momentum flux 
ratio which results froa using the average velocity exiting froa the 
sharp edge orifice. The velocity profile at the exit of the nozzle 
used by Kaaotani and Greber (1971) is iiCaerently acre uniform than 
that at the exit of a sharp edge orifice and the momentum flux ratio, 
based on an average velocity, will be less for their nozzle. Another 
important difference between the two experiments is the cross stream 
boundary layer thickness compared with the jet exit diameter; 5/d < .1 
for Kamotani and Greber, 5/d * 1,3 for the present experiment. 

3.5 Comparison of undisturbed with disturbed cross flow conditions 


Figures 28-33 are contours of mean isotherms for the case where 
the jet exits into a disturbed cross flow. The shapes of these con- 
tours are very similar to those for the undisturbed cross flow case, 
(Figures 22-27); however, the maximum mean values for the disturbed 
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cross flow cate are significantly lower than those for the undisturbed 
cross flow case. Additional differences are reflected in the contours 
of constant tomperature standard deviations (Figures 44-55) and the 
associated histograms (Figures 35-42). 

The contours of constant standard deviation show dramatic differ- 
ences between the disturbed and undisturbed cross flow conditions. 
The values of the standard deviations are larger for the disturbed 
cross flow condition and the temperature standard deviation field is 
larger than that for the undisturbed case. 

The histograms reveal a great deal of information regarding the 
differences between the undisturbed and disturbed cross flow cases. 
For example, compare figures 35 and 39, These two figures represent 
the histograms for thermocouples located near the maximum contour 
level for tho specified J value and for the undisturbed and disturbed 
cross flow conditions respectively. The bimodal characteristic for 
the disturbed cross flow case, shows that ambient (cross stream) tem- 
perature fluid occupied this point in space for a significant fraction 
of time. Similar histograms exist for the larger J values; see Figure 
37 and 41. These suggest that tho jet trajectory excursions are a 
significant fraction of the jet width for the disturbed cross flow 
case&. 

Six instantaneous time series, each approximately .8 seconds in 
length, were taken and stored for the Tj-x^dl ,1*C overheat cases: 
J=21 and 84, disturbed and undisturbed cross flow conditions. These 
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ITS. 



J 


®t 

c.g. 

K.t 

undisturbed e.f. 

21 

.341 

.233 

.067 

undisturbed c.f. 

84 

.233 

.233 

-0.13 

disturbed c.f. 

21 

1.84 

1.43 

-0.33 

disturbed c.f. 

84 

1.16 

1.12 

-0,43 


TABLE II. 2 Cross stresa disturbance effect on centroid of 
temperature field 


time series were formed from periodic (At*0 , 64msec ) and simultaneous 
samples of the 64 thermocouples. The centroid (s c g <t c g ) of the 
instantaneous temperature field was determined using equation II. 7 and 
representative time series for the four cases are shown in Figures 
56-59. (Note that the two coordinates of the centroid: * c ,g. and 

^c.g.* tre shown separately). The increased movement of the jet 
column for the disturbed cross flow case is clearly seen in these fig- 
ures. The following table shows the results using all six of the time 
series for the disturbed cross flow case and three of the time series 
for the undisturbed cross flow case. 


The correlation coefficient (K st =it/ (or s a t ) results presented in 
the above table show that the Y and Z motions of the jet are signifi- 
cantly correlated for the disturbed cross flow case. This correlation 
is considered to be representative of the cross flow itself rather the 
dynamics of a jet in a disturbed cross flow. That is, the correla- 
tion, that is suggested by the K st values, would be expected from the 
large-scale secondary motion, for the related problem of a bounded 
jet, as documented by Holdeman and Foss (19). The time series 
records, Figures 56-59, indicate a low frequency translation of the 
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centroid. This it qualitatively similar to the (hot-wire anemometer 
derived) energy spectrum results presented in Figure 16. 

Instantaneous isotherms can also be generated from the time ser- 
ies taken at the high sampling rate. A selected few of these 
isotherms are shown in Figures 60-67, These isotherms were selected 
based on the centroid time series shown in Figures 56-59 . Referring 
to figures 36 and 57 (J*21), one can see that at approximately .2 and 
.65 seconds there are peak excursions of the jet centerline for the 
disturbed cross flow case. The isotherms for both the disturbed and 
undisturbed cross flow case at these two instants in time are shown in 
Figures 60-63. At the discrete points in time of approximately .3 and 
.38 seconds, figures 58 and 59 <J»84) show peak excursions of the jet 
centerline. The corresponding contours can be seen in Figures 64-67. 

The migrating jet column is clearly seen when comparing the dis- 
turbed cross flow isotherms. For the disturbed cross flow, J*84, 
contours shown in Figures 66 and 67, the excursion is so great that 
the resulting two temperature fields hardly overlap. This is also 
true for the J=21 contours seen in Figures 62 and 63. The migration 
of the jet column for the undisturbed cross flow case is not as severe 
as can be seen in Figures 60, 61. 64. and 65. 

Comparing the instantaneous contours with the mean contours, one 
can see th?t the differences are greater for the disturbed cross flow 
case. The maximum values for the instantaneous contours of the dis- 
turbed cross flow case are much larger than those of the mean 
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contours. The maxinus instantaneous values for the disturbed cross 
flow case do uot differ significantly fro* those for the undisturbed 
cross flow case. The mean temperature values are such lower however, 
because of the jet centerline excursions. 
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CHAPTER 4 


CONCLUSIONS 


The following conclusions are based on the results of the present 
s tudy . 


i) Substantial diffusion of thermal energy has occurred between 
the exit plane and measurement location, since the maximum tempera- 
tures recorded were <.25 (Tj et -T m ) for both the undisturbed and 
disturbed cross flow cases. 

ii) The present data allow the sensitivity of the jet penetra- 
tion, to the magnitude of the momentum flux, to be quantified. In 
general, the disturbed cross stream case is less sensitive to the 
momentum flux ratioalbeit the differences, between these two condi- 
tions, are neither large nor simply expressed. 

iii) The present results show that the centroid of a heated jet 
has much larger excursions, for the jet exiting into a highly dis- 
turbed cross flow than for a jet exiting into an undisturbed cross 
flow. 
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iv) For both the undisturbed tod disturoed cron flow cases. an 
inatantaneous isotherm pattern can be significantly different from the 
mean isotherm pattern and this difference is more dramatic for the 
disturbed cross flow case. Specif icallyi the instantaneous maximum 
temperatures for the disturbed flow are similar in magnitude to those 
of the undisturbed case whereas the maximum mean temperature is sub- 
stantially less for the disturbed cross flow compared with the 
undisturbed cross flow case. 
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b) Detail of Temperature Sensing Junction 
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Fv'jure 3. 1 of 18 Thermocouple Modular Elements 
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Figure 5. Schematic Representation of One of Four Electronic 
Processing Boards 














Figure 6. Detail of a Sample and Hold Circuit 
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Figure 7. Experimental Configuration for the Time Constant Evaluation 



Figure 8. Schematic Representation of a Thermocouples Response to the 
Impulsive Removal of a Thermal Input 
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Figure 11. Schematic Representation of Jet Trajectory, Showing Lab 
Coordinates (x, y, z) and Jet Coordinates ( S, n. ;) 




Figure 12. * Schema tic Representation of : the TSFL Wind Tunnel 






Figure 13. Normalized Velocity Distribution of Shear Layer at 
Jet Exit Location 
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Figure 19. Placement of Thermocouple Array in Flow Facility 




Figure 20. Cross Stream Boundary Layer for Undisturbed Cross 
Flow Case 
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Figure 22. Undisturbed Cross Flow (u.c.f. ) Mean Temperature 
Isotherms j'= 17 » Tj et - T oj = 22.2°C 

T ' T CC 

NOTE: Contours shown are f- X 1000 s/d = 14.1 (+ .7) 

jet - =° 

t/d = 4.3(+ .5) - y/d 









Figure 24. u.c.f. Mean Temperature Isotherms 
NOTE: s/d = 14.1{+ .7) - z/d, t/d = 4.3(+ .S 
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Figure 27. u.c.f „ Mean Temperature Isotherms: J = 84, T. gt - T m = 61.1°C 


NOTE: s/d = 17.9(+ .7) - z/d, t/d = 4.3 (+ .5) - y/d 
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Figure 30. d.c.f. Mean Temperature Isotherms: J =21, T 

NOTE: s/d = 14.K+ .7) - z/d, t/d = 4.2(+ .5) - y/d 
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Figure 32. d.c.f. Mean Temperature Isotherms: 3 = 1 
NOTE: s/d = 16.8(+ .7) - z/d, t./d = 4.2(+ .5) - y/d 
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Figure 34 Thermal Centroids as a Function of the Specific 
Momentum Flux Ratio 

Note: Solid Curve added for visual reference only. 
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Figure 37. Temperature Histogram for 
ThermocouDle at s/d = 4.7 
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Figure 41. Temperature Histogram for d.c.f 
Thermocouple at s/d = 3.87, t/d 
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Figure 45 


u.c.f, Contours of Constant Temperature Standard 
Deviation: T. t = 41.7°C 


NOTE: Values shown are T OOOX[T_{ s ,t) - T 1 / (T . - T ) 

u 00 rms j °° 

s/d and t/d as defined for the mean temperature isotherms 





Figure 46. u.c.f. Contours of Constant Temperature Standard 
Deviation: J=21, T^ et - - 61.1°C 

NOTE: Values shown are 1000X[T c (s,t) - Tj^/tT^ - Tj 

s/d and t/d as defined for the mean temperature isotherms 






Figure 47. u.c.f. Contours of Constant Temperat: 

Deviation: j=67, TV* - T = 22.2' 

jet =° 

NOTE: Values shown are 1000X[T„(s ,t) - T 1 /(' 

*■ c' •=°-‘ rms 

s/d and t/d as defined for the mean tempe 
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Figure 48. u.c.f. Contours of Constant Temperature Standard 
Deviation: J=76, Tj et - T m = 41.7°C 

NOTE: Values shown are 1000X[T c (s,t) - ■T ro 3 rms /(Tj - Tj 

s/d and t/d as defined for the mean temperature isotherms 






Figure 50. Disturbed Cross Flow (d.c.f.) Contours of Constant 

Temperature Standard Deviation: j = 17 T . . - T = 22.2°C 

•jet » 

NOTE: Values shown are T000X[T c (s ,t) - T co ] rrns /(Tj - Tj 

s/d and t/d as defined for the mean temperature isotherms 
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Figure 51. d.c.f. Contours of Constant Temperature Standard 

Deviation: J = 19, T- . - T = 41.7°C 

jet ® 


NOTE: Values shown are 1 000X[T c ( s ,t) - Tj rms /1\ - Tj 

s/d and t/d as defined for the mean temperature isotherms 




Figure 52. d.c.f. Contours of Constant Temperature Standard 
Deviation: J = 21,T- et - = 61.1°C 

NOTE: Values shown are 100QX[T e (s,t) - 'f m 2 rms /T - - Tj 

s/d and t/d as defined for the mean temperature isotherms 




Figure 53. d.c.f. Contour of Constant Temperature Standard 
Deviation: J = 67, Tj et • T a s 22.2°C 


NOTE: Values shown are 1000X[T c (s ,t) - T oD l rms /T j . - Tj 

s/d and t/d as defined for the mean temperature isotherms 
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Figure 56. Variation of s Component of Centroid with Time 






Figure 58. Variation of s Component of Centroid with Time 
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Figure 61. u.c.f. Instantaneous Temperature Isotherms 

J = 21 T. 4 . - T = 61.1°C, Time = .65 sec 
’ jet ” 


NOTE: Contours are j - 2 - — X 1000, s/d and t/d as defined for 

'jet - 

the mean temperature isotherms 









Figure 63. d.c.f. Instantaneous Temperature Isotherms 
0 = 21, Tj e _j. - ® 61.1°C, Time = .65 sec 

V T » 

MOTE*. Contours are j- j- X 1000, s/d and t/d as defined for 

‘jet ~ 

the mean temperature isotherms 





Figure 65. u.c.f. Instantaneous Temperature Isotherms 
J * 84, Tj et - T w = 61.1°C, Time = .58 sec. 


NOTE: Contours are T 0 X 1000, s/d and t/d as defined for 

jet » 

the mean temperature isotherms 




Figure 66. Disturbed Cross Flow (d.c.f.) Instantaneous Temperature 
Isotherms, J = 84, T - T w * 61.1°C, Time = .3 sec. 


T -T 

0 03 


NOTE: Contours are f ■ _ T X 1000, s/d and t/d as defined for 

jet - « 

the mean temperature isotherms 




Figure 67. d.c.f. Instantaneous Temperature Isotherms 
J = 84, Tj et - T m = 61 .1°C, Time = .58 sec. 


NOTE: Contours are ■ X 1000, s/d and t/d as defined for 

jet » 

the mean temperature isotherms 
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